Pseudogap, competing order and coexistence of staggered flux and d-wave pairing in 

high-temperature superconductors 
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We study the t-J-VmoAe\ of a doped Mott insulator in connection to high-Tc superconductors. 
The nearest neighbor Coulomb interaction (V) is treated quantum mechanically on equal footing as 
the antiferromagnetic exchange interaction (J). Motivated by the 5*17(2) symmetry at half-filling, we 
construct a large-N theory which allows a systematic study of the interplay between staggered fiux 
order and superconductivity upon doping. We solve the model in the large-N limit and obtain the 
ground state properties and the phase diagram as a function of doping. We discuss the competition 
and the coexistence of the staggered flux and the d-wave superconductivity in the underdoped regime 
and the disappearance of superconductivity in the overdoped regime. 
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The physics of doped Mott insulators has played a cen- 
tral role in the study of high-Tc superconductors 0, 01 ■ 
Traditionally, the Mott insulating state is derived from 
the strong onsite Coulomb repulsion and naturally ac- 
counts for the antiferromagnetic (AF) ordered state at 
half- filling 1]. There is increasing evidence that, due to 
large quantum fluctuations in spin-1/2 systems, several 
exotic quantum spin liquid states are likely hidden under 
the AF Mott insulating ground state. These states are 
competitive in energy and form a basis for the resonating 
valence bond route to superconductivity 0| . Several im- 
portant, yet unresolved problems remain: Is the d-wave 
superconductor (dSC) the ground state when the long- 
range AF order is destroyed by the coherent motion of 
the doped holes? Does the competition between the dSC 
and other forms of order account for the pseudogap phe- 
nomenon in the underdoped regime? Does the physics 
of doped Mott insulator explain the disappearance of su- 
perconductivity in the overdoped regime? 

In this paper we demonstrate that long-range Coulomb 
interactions, not crucial for discussing the Mott in- 
sulating state, play an important role in addressing 
these issues associated with doped Mott insulators. Be- 
cause charge fluctuations are limited by the density of 
doped holes, long-range Coulomb interactions have been 
thought to be not important near the stoichiometric Mott 
limit. However, since the strength of the interaction is 
large compared to the Fermi energy, its effects become 
quickly significant with doping. We consider the mi- 
croscopic t-J-VxnodeX where V stands for the nearest- 
neighbor (NN) Coulomb repulsion. As in the t-Jmodel, 
the constraint of no double occupancy leads to an in- 
sulating state at half-flUing. It is well known that the 
d-wave pairing state and the staggered flux phase (SFP) 
0,0,0, identical at half-filhng owing to an SU{2) gauge 
symmetry 0, are two competing phases at low doping. 
The SFP has a plaquette flux density wave (orbital cur- 
rent) order with d-wave symmetry. It is a metallic state 
with hole pockets and a pseudogap of the same symmetry 



as the dSC. Thus the SFP serves as a candidate pseu- 
dogap phase for underdoped cuprates, much as in the 
d-density wave scenario 'J']. However, in the standard t- 
Jmodel, the ground state becomes the dSC upon doping 
0|- The SFP was found to be unstable to superconduc- 
tivity [lol ITU and only exists at finite temperatures in a 
very narrow regime near half-filling [l^ . 

We show that the NN Coulomb interaction introduces 
competing and coexisting order and changes significantly 
the structure of the phase diagram. It is crucial to treat 
the V^-term quantum mechanically, on equal footing as 
the exchange J-term, since both pairing and screening 
originate at the same energy scale. Motivated by the 
SU{2) symmetry at half-filling, we construct an SU{N, q) 
{q =quaternion) generalization of the t-J-VmoAel. While 
the SU{2N) theory |l3l favors the normal state and 
the SP{2N) theory '14'| favors the superconducting (SC) 
state in the large-N limit, the SU {N, q) theory accommo- 
dates both staggered flux (SF) and SC order, and allows 
a systematic study of their interplay. In the large-N limit, 
we find that the SF order coexists with the dSC below a 
critical doping Xc where a quantum phase transition to a 
pure dSC phase takes place. The SC phase is destroyed in 
the overdoped regime beyond Xp > Xc- At finite temper- 
atures, the onset of the SF order is at Tsf > 7c, giving 
rise to a pseudogap SFP above the SC transition temper- 
ature Tc- We discuss the phase diagram, the signatures 
of the coexistence and the pseudogap phase. 

We start with the large-N generalization of the t-J- 
l^model on a square lattice, H = Htj + Hy, 
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The operator c|g. creates an electron with spin cr = 
1,1; ■ ■ ■ ; N, N at site i under the constraint of no dou- 
ble occupancy for each spin flavor, rii — cl^Cia- < N. The 
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NN Coulomb repulsion between the electrons is given by 
Vc-, which is grouped together with — {J/A)ninj from the 
exchange term and represented by the V-ievva in Eq. |(2Jl. 
The usual t-J-Fmodel corresponds to iV = 1. The pro- 
jected Hilbert space can be treated in the slave boson 
formalism by writing — fj^hi, where //^ is a spin- 
carrying fermion and bi a spinless boson keeping track of 



empty sites, with the constraint fj^fia + h\h.. 
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Let us introduce the pseudo-spinors, = (V'iTjV'ii) 
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and V'zj, — Ci/'it where C = iay®\. Here Iq,^ is the iden- 
tity matrix in the flavor space and {<yy)TT' is the Pauli 
matrix acting on the spins. Under a unitary transfor- 
mation gi, ipi-^ gi'ipii and ipj^ ''^lidl^ provided that 
glCg^ = C, i.e. g, e SU{N,q). The SU{N,q) fermion 
bilinear Uij — ^i^*] transforms as Uij gtUijg^y To 
construct the large-N theory, we form the flavor singlet 
by tracing over the flavor indices. 



where g^ = J/4: + XiXjV and gA = J/^ — XiXjV are the 
fermion coupling constants to the bond and paring order 
parameters respectively. The effect of V is to reduce 
pairing while enhance the SF ordering. We focus on the 
uniform solutions with Xi = x, iXi — A, Xij = Xi + 
i(_l)«.+jy^2 and Ay ^ As + i(-l)^«+J«Ad where xi,2 
and Ag.d are real. The flux per plaquette is given by 
(j) — 4tan~^(x2/xi)- and Ad are the s-wave and 
d-wave pairing order parameters, respectively. The free 
energy per site is obtained analytically, 

- ^TN^ ln[2cosh(£;fes/2T)] 



k,s=± 
|2 I „ U,|2n 
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where fj, = iif ^ X. The fc-sum is over the reduced zone 
due to the doubling of the unit cell, and E^s describe the 
two branches of the quasiparticle (QP) dispersion 
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which is a quaternion with elements: Xij = X^qI/jL/jq + 
fLfjs) and Ay = EaifLfja - flafjJ- They represent 
the bond (SF) and pairing degrees of freedom. In terms 
of u, the Hamiltonian can be written as. 



Hu = - —T 
2N ^ 



hb] V-t (1 + <^z)i'jr + h.c. 



(5) 
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Y b] bi b] bj tr (u^ CT^ Uij ) . 



(6) 



Here tr stands for the trace in spinor space, \i is the La- 
grange multiplier enforcing the local constraints, and /i/ 
fixes the average fermion density. Note that for general 
iV, Hy cannot be written solely in terms of four fermions 
using the constraints. 

In the large-N limit, the bosons condense at low tem- 
peratures (&iL= ^/Nxi, where Xi is the local doping con- 
centration [13. Note that since charge fluctuations are 
suppressed, the V-term vanishes with vanishing boson 
amplitude on approaching half-filling where the SU{2) 
symmetry is restored. The interactions are decoupled by 
introducing the saddle point values Xij — iXij)/^ and 
Aij = {Aij)/N. The saddle-point Hamiltonian becomes 



H^j = - ty/^Xij - 9xx\jXi] - 5aAJ^- Ay + h.c. 

+ iV(gJxyf +5A|Ayf), (7) 



with 7± = coskx ± cosfcy, — 2{xt + gxXi)l+^ and 
[5k — 2g^X2l-- Minimizing the free energy in Eq. ((SJ, 
we obtain the self-consistent equations. 
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In the absence of the NN Coulomb interaction, V = and 
9x — 9^- The equations for X2 and A^ imply that they 
cannot be simultaneously nonzero except at half-filling 
(m = 0) HO] ■ The NN Coulomb interaction > leads 
to g^ > away from half-filling, making it possible for 
the coexistence of the two types of order. 

We solve the self-consistent equations in (O using J 
as the energy unit and t — 2J, V — 2J. The s-wave 
pairing As is found to be zero. In Fig. 1, A^, xi and 
X2 are plotted at T = as a function of doping x. The 
results for = is also plotted for comparison. The NN 
Coulomb interaction introduces two important features, 
(i) The dSC order develops upon doping, but is destroyed 
beyond a critical doping Xp ~ 0.26, which is smaller than 
the naive estimate from the condition ^a = 0. The sup- 
pression of pairing by V in the overdoped regime is con- 
sistent with recent quantum Monte Carlo studies of the 
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FIG. 1: The order parameters versus doping at T = 0. 

t-U-V model (ii) The SF order develops below a crit- 
ical doping Xc — 0.1 and coexists with the dSC. There 
is a reduction of concomitant with the development 
of the SF (X2) below Xc, indicative of competing order. 
The values of Xc and Xp depend on the band structure. A 
next NN hopping t' <Q moves both to higher values. Al- 
though the d-wave pairing amplitude A^j extrapolates to 
a nonzero value at half-filling, the two states are identical 
by the SU{2) gauge symmetry. 

In Fig. 2, the phase diagram of the t-J-Vmode\ is 
shown in the plane of temperature T versus doping x. 
The phase boundaries are the temperatures Tsf and Tc 
at which the SF and dSC order parameters vanish, re- 
spectively. The topology of the phase diagram resembles 
the one proposed for the cuprate high- Tc superconductors 
by invoking competing order and quantum criticality 
The dSC sets in below Xp. The transition temperature Tc 
rises with decreasing doping and reaches its optimal value 
somewhat below the quantum critical point (QCP) at Xc 
where coexisting SF order emerges. The Tsf line shows 
a bulge when crossing the Tc line at x'^ ~ 0.12, making 
it possible for a thermal transition from the dSC to the 
coexistence phase for Xc < x < x'^. The phase bound- 
ary rises almost vertically at low temperatures, and the 
transition is weakly first order as is evidenced by the 
small discontinuity in the order parameters (Fig. 1). For 
X < x'^, Tc < Tsf, giving rise to a pseudogap SFP with 
orbital currents for Tc < T < Tsf in the underdoped 
regime. The phase diagram and the nature of the pseu- 
dogap phase are different from those in the SU{2) for- 
mulation of the i- Jmodel % In that theory, the SU{2) 
symmetry at half-filling is extended to finite doping via 
the introduction of a slave boson doublet. The pseudo- 
gap phase is a fluctuating state between gauge equivalent 
SFP and the dSC. The static SF order develops only in- 
side vortices where SC order is suppressed by the external 
magnetic field 17]. The present theory explores the con- 
sequences of symmetry breaking by doping, through the 
NN Coulomb interaction in particular, which stabilizes 
the static order of the SF in the pseudogap phase. 

To illustrate the electronic structures, we calculate the 
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FIG. 2: The large-N phase diagram of the t-J-Vmode\ show- 
ing a QCP at Xc, a pure dSC, a coexistence phase of dSC 
and SF, and a finite temperature pseudogap SF phase. Inset: 
T = single-particle gap Aq, condensation energy Ai?c(x50), 
and specific heat coefficient jump Ajc versus doping. 

tunneling density of states (LDOS) measured by STM 

/•OO 

iVM = -Re/ dte'-^'0{t){{cl{O),c,„{t)}). (10) 
Jo 

In Fig. 3, we show the typical LDOS and the QP disper- 
sion. Fig. 3(a) corresponds to the pure dSC at x = 0.14 
and T = 0. The LDOS shows d-wave gap features. The 
two peaks near the gap edges come from the large state 
density at the points marked as X at (tt, 0) and P at the 
nearby minimum on the dispersion, and their symmetry 
equivalent points. Note that the distance between nei- 
ther set of peaks equals to 2 A where A = A^ is the 
maximum of the d-wave gap function. One must thus be 
careful with extracting the SC gap from the peak-to-peak 
distance even with high STM resolution. The LDOS in 
the coexistence phase is shown in Fig. 3(c) at 2: = 0.06 
and T = 0. Compared to Fig. 3(a) in the dSC, there is 
an emergent in-gap resonance peak (marked A) particle- 
hole reflected on both sides of zero-bias. It arises from 
the resonant scattering of the QP from state k to state 
k + (tt, tt) due to period doubling. A similar situation 
arises in the SF state around a superconducting vortex 
in the SU(2) formulation of the t-J model [ll|. As shown 
in Fig. 3(d) for the s = + branch, at the resonant energy, 
the distorted ellipses touch at point A, the low-energy 
minimum in the dispersion plotted in Fig. 3(c). These 
resonance peaks, albeit small in weight, are ubiquitous 
signatures of coexisting order and should be detectable 
by high resolution STM. Near the gap edges, there are 
four sets of peaks (two are prominent) due to the dou- 
bling of the unit cell. They are labeled by X± and P± 
and marked on the QP dispersion. It is instructive to 
compare these to the LDOS in the pseudogap SFP at 
X = 0.06 and T = 0.08J shown in Fig. 3(b). The low en- 
ergy resonance peaks disappear and the LDOS shifts to 
the unoccupied side. Interestingly, the minimum of the 
LDOS in the SFP coincides with the M point next to the 
in-gap resonance in the coexistence phase. The resonance 
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FIG. 3: The LDOS and the QP dispersion in dSC (a), SF 
phase (b), coexistence phase of SF and dSC (c). (d) Equipo- 
tential contours at = 0.03, 0.04, 0.043, 0.05J, showing reso- 
nant scattering at uja ~ 0.043J. 

energy and the spectral shift depend on the band struc- 
ture and doping and are on the order of 0.043 J ~ GmcV 
for the parameters used in Fig. 3. 

An important consequence of the competing order is 
the emergence of two energy gaps which can be, how- 
ever, difficult to observe at low temperatures [19j- In 
the inset of Fig. 2, the gap to single-particle excitations 
Aq obtained from the LDOS at T = is plotted versus 
doping, which qualitatively agrees with the low temper- 
ature ARPES data. It interpolates between a geometric 
average of the gaps associated with the SF and SC or- 
der for X < Xc and the gap of the pure dSC for x > Xc 
and closes at Xp. The only sign of competing order in 
Aq is the mild cusp-like singularity near Xc- Signatures 
of competing order are more apparent in thermodynamic 
quantities same figure, the jump Aj^ 

in the specific heat coefficient 7(T) = Cp{T)/T at 
and the condensation energy per site AEc at T = are 
shown as a function of doping. Ajc rises with decreas- 
ing doping on the overdoped side. With the opening of 
the pseudogap due to SF order, the entropy in the nor- 
mal state is reduced, Ajc turns around at x'^ and drops 
continuously to small values in the underdoped regime, 
in q ualitative agreement with specific heat measurements 
[l9j . The condensation energy AEc grows almost linearly 
with doping in the underdoped regime and peaks around 
Xc- The closeness in the free energy between the SFP 
and dSC, the smallness of the condensation energy, to- 
gether with a small superfluid density that scales with x 
promote the appearance of cheap vortices in the under- 



doped regime as emphasized in the SU{2) theory [l^ . 
except that static orbital currents in the present theory 
live in the bulk as well as inside the vortex. 

In summary, we have presented a microscopic large-N 
theory for the competition between the exchange and the 
finite-range Coulomb interactions in doped Mott insula- 
tors. The i-J-V^model shows competing and coexisting 
order of SF and dSC which can be viewed as a product 
of SU{2) symmetry breaking by doping through the NN 
Coulomb interaction. The phase diagram exhibits a QCP 
separating the coexistence phase from a pure dSC, a finite 
temperature pseudogap SFP in the underdoped regime, 
and complete suppression of superconductivity beyond 
an overdoping. At very low doping, the phase diagram 
does not match that of the cuprates. The discrepancy 
is most obvious in the nonzero Tc that extends to in- 
finitesimal doping, whereas the SC state in the cuprates 
disappears at a finite doping below which a spin glass 
phase prevails. The lack of destruction of the dSC at low 
doping is because the renormalized NN Coulomb interac- 
tion must scale to zero at half-filling due to the constraint 
forbidding charge fluctuations. We expect the effects of 
doping induced inhomogeneity [ij due to poor screening 
of the dopant ionic potential to become very important in 
this regime As a result, the local doping concentra- 
tion (LDC) fluctuates dramatically and its interplay with 
AF may lead to the destruction of the SC state and the 
emergence of spin glass and insulating behaviors. Fur- 
thermore, the onset of SF order will fluctuate spatially 
with the LDC, smearing out the sharp finite temperature 
transition to the pseudogap SFP. 
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